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AN INTRODUCTION TO ENGINEERING MATHEMATICS, 
WITH REFERENCE TO HEAVISIDE’S 
OPERATIONAL CALCULUS. 


(With Examples of an Electrical Nature) 
BY 


Ronatp J. Brrxrnsuaw, A.R.T.C.S., Grap.LE.E., Grav.I.Ec.E., 
M.J.Inst.E. 


INTRODUCTION. 


The following section is merely a brief outline of the various 
methods used in the calculus. 

For reasons of space-economy, the subject is not treated at 
great length, as many excellent text-books on this branch of 
mathematics, already exist. 

A few of these books are mentioned in the reading list at the 
end of this pamphlet. 


BRIEF RESUME OF THE CALCULUS. 
Differentiation of y =f (x). 


Let y =x 
Imagine y given a small increase Sy, and x a small increase x. 
Thus y + dy = (x + 8x)8 
= 8+ 3x2 6x+3x (dx)? + (dx)8 
Sy = «3+ 3x? Sx + 3x (5x)? + (8x)8—x? (=) 


3x? 8a + 3x (dx)? + (Sx)3 
As 8x is very small, then powers of this quantity may be 
neglected. 


dy 
= Oa 
3x 3x 
In the limit as 4 -—> 0 
dy d (fr) : 
a = 7E = 3% 


In general if y =a, es = nxe-t 
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Introducing a constant 6, where 6 is not a function of x, then if 
y = bx, 2 = nbz 
The term dy/dx is the first differential coefficient, or first de- 
rivative, and sometimes denoted by f« (x). 


Successive differentiation is carried out in exactly the same 
manner. 


Thus y = f(x) =a" 
Bape) = nem 
= = f™ (x) = 2 (w—1) xe? 
Xp (x) = 0 (=I) (9-2) 204 
oa =f" (@) =n (n—1) (n—2) (n—3) 40-4 


“ =f" (x) =» (n—1) (n—2) (u—3) (n—4) xn-5 


etc. 
Suppose y = #", where w is a function of x. This is known 
as a function of a function. 


Then ay ay aE 


RATES OF CHANGE. 


If y is a function of time, ¢ (=f ()), then the value of yata 
time ¢ + 8f is y + Sy. 


During the time interval from ¢ to ¢+8¢, y changes at the rate 


of dy/8t. im 
As &¢ tends to the limit, & = 0 
by dy 
e becomes oF 


Referring to Fig. 1, let Q be a point on the line OX, moving in 
a direction O to X, then at a time ¢, the velocity of Q = dx/dt. 
Thus V = dx/dt. ah 
The acceleration is the rate of change of velocity. With 
respect to time = du/dt 
dv ax 


Acceleration = a = og 


. 
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) ET sie 
Fig. 1.—To illustrate rate of change. 


The following is a list of some important differential coefficients. 
y = e* [e is the base of natural logarithms]. 


eet ee oat Gh zm 
e= +X + [2 ? [30% (4 ae oe [n 
ane Pr x3 * gn 
= (2) =10 + 1 Set Re we ot gy ee 
dx (e) for et |2 a |3 ‘ |n 

d " dy 

AM ex\ i ee Gk 

aa piace 


y =v [y, v and u being functions of x]. 
Let each function be given a small increment. ; 


Then y + Ay = (u + Au) (v + Av) * 
y+ Dy = wo + uAv + vAu + Au Av 
y = uv 
Ay u Av vAu Au Av 
= + aS 
Ax Ax Ax Ax 
MUHA 
as Ot as Ax .—> 0 
dy Lyra du 
dx Ta ae d. 
hoa 
y ae Si {obeying the same premise as before]. 
v : : ; 
u+ At % 
pia ae y . Rta 
f w+ Au u- 
49 a) + Av Ders 
_ uv + vAu — uw — uv 
ir v (v + Lo) 
ie vAu — uv 
a (v + Av) 


6 AN INTRODUCTION TO ENGINEERING MATHEMATICS 


Nu Lo 
Ay a Ax iss Ax 
Ax v (v + Av) 
dy Ay . z - 

’ = > 0 
dx _ Le $3) O48 
,, du dv ; 

dy Te dx 
dx v2 


PARTIAL DIFFERENTIAL. 
It is found in various applications of mathematics, that functions 
of more than one variable occur frequently.’; 
. Suppose « =f (x.y) 
In this type of differentiation a different differential coefficient , 
symbol is used ; it is a “curly” d, thus :—é. 
To find f* (¥4), each variable is differentiated separately. 
Let « = ax® + bumy + CxyP 
Ou 
pe axl + mbxmly + CypP 
3 = bx + pCxyrt 


INTEGRATION. 


Integration is the reverse of differentiation. Let f(x) be a 
function of x. 


> fe -fo- [rea 


f (x) is the integrand and x is the variable of integration. This 
type of integral is known as indefinite. The insertion of limits 
transforms it into a definite integral, e.g., 

. x=a 
| FT (*) dx 
xoO 


If b is a constant (z.e., not a function of x), then’ . 


[ bf (x) dx = | rt) de 
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Sal 


y c ype 
Example: -} x°dx = 


[ provided n is not = -1| 


SUCCESSIVE INTEGRATION. - 


Suppose the density # of a plate varies as its dimensions x and y. 
Then # may be expressed as a function of these two quantities, 
b =f (xy-) : 


The evaluation of this quantity leads to a double integral 


[. i f (wy) dx dy 


Similarly, if the density varies also as the thickness ¢; then a 
triple integral results. 


thus; ;b -d cf ° 
| | f (xy 2) dx dy dt 
AEN Ce 
: . Successive integration may also be applied to a, particular 
function. 
In the case of beam problems, it is known that : . 
CURVATURE = —— 
TURE El I 
SLOPE etl are 
nl nee 
DEFLECTION = | |= M dx dx 
EI 


(Problems of a similar nature will be discussed in. more detail ' 
in another part of this section). 


INTEGRATION BY PARTS. 
It has been shown that : 


d. dy du 4 
Pra A ta Er 2) 
| : —4 wy =" Be — a uv ¥ 
ine - dx ee ay dx ‘ 
Multiplying by —1 and integrating, 


dv du 
[ute wv — | a mn” 
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DIFFERENTIAL EQUATIONS. 


These form a very important part, in the solving of many pro- 
blems arising in electricity, mechanics, chemistry, physics, etc. 
An equation containing derivatives of a variable is called a differén- 
tial equation. . 

The order of a differential equation depends upon the power 
to which the highest derivative is raised. 

The degree of a differential equation is the degree of the highest 
differential coefficient. 

As an illustration of these definitions, consider the following 
equations :-— : 

dy d*y dy 


L b + = 
be eiasr pia ieee ae FA 


This is of the third order, and first degree. 


dy . dy 
b > = 
) ie +m ae f (x) 
This is of the second order, and first degree. 
2 Berns ny) : 
a —_— t — 
¢) ap thf 


This is of the first anele, and ee degree. 
a’y dy ; 
d (2) BEA ta no 
DUN rs) i Paes § oS) 
This is of the second order, and second degree. 
dy) @y ty 
iy (¢) (aa) Mae ae nals 
This is of-the fourth order, and third degree. 


All the above equations are ordinary differential equations. rAS 
- partial differential equation is one involving partial differential 


coefficients. 


a a ay 
oy + BP sia 22 i 
em? om Ox? Ox 


Cheb CH 


‘Types of Differential Equations. ° 
(i) Exact DIFFERENTIAL EQUATIONS. 


Integrating, x.y = C (where C is a constant). 


, 
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Example — : : 
1, 5 +Ri=E 


Multiplying by eR/!-* 


ee oRlb-t WRG neki ails are 
a eRIL.t oy Bi. enn = Eek. 
yi 2 $ | eRIL-t 
Integrating, a. eRL-t = 7 e-RIL.t 4 K 
Multiplying by eRIL-t 
2 e + K.e-RL-t 


(#) SOLUTION By SEPARATION OF VARIABLES. 


ey BINS 
1g a +°Ri = E 
di: E-Rt 
CB Ty ON 
Pe 4 Fel 
“peg ane = - . dt = 
ae als 1a°4 
1 R 
mm SS EB 
: “E/R+i' 13 
4 E 
Integrating, log [- mt i] =- —.t+A 
" e. E +4 = er eRIL-t 
R 
let e* “>= K 
i= Bz + Ke-RL.t 


—- «2 
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(#11) Homocenrous Equations. 


gant) 
let vad y = ox 
mote, & 
vt. Safe 
do 


iy = —.dx 
i OeD % 
Integrating 
{ 1g - dx = log x + C 
Palaa y= |. de = log x + 
Example. 
dy x+y 
dx a 2 xy 
y 
(Ci aaa: 
ae 
dy | qd 1+ 
dik t aE Qu 
dv Lev? 1 -v? 
i ee ne ea hes 
dx 20 B 20 
20° 1 
I, dy =—., 
ieviet eas x fp 
rath 6 - log (1-v’) = log x + log C 
. (l-v) +¥ +C =0 
ee — y+ x3, 
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(iv) LryEar Equations. 


dy 
pas Py = 
Pep ae Q 
Reducing this equation by equating Q to zero, we have : 
dy 
= Py = 0 
CE ae 
oy = — Pivdx 
Integrating, log. y = -— JP.dx + log. C 
i) y = efPdeC 
€= y efP.dx 


Differentiating, 
dy * 
sP.dx | = 0 
2 Ee tare ] 


also, of? ae + Py | = Qe/P-d 


d 


= (y-efP a) = Qe/P-a 
Integrating, ' 
Hyer uede . = | QefP.de dx + i 
y = Peds [: [ Qer- w+ K | 
Example. ‘ x 


Consider the equation, L di/dt + Ri = E, the solution of which 


has been found previously : 
E . 
= eaRibet [le GER/Let dt as K| 
E is constant (D.C. case). 


If at ¢=0, 2=0, 
then : 


ee a - Ret | 


il 
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(v) Lingzar DirFerentraL EQuaTion HAVING CONSTANT 


COEFFICIENTS. . 
aty dy 
a. a? Obl Peter eda 


* The following may be noted in connection with the above | 
equation : 


(a) The coefficients a, b and ¢ are constant. 

(6) The equation is of the second order. ; 

(c) The equation is linear, i.e., the highest derivative is of the 
first degree. i , 


Example. 
ay dy 
Rae ce oe 1 = 
FE ede ais 1 
Let y be of the form e™ A | 
dy mt ay 2 gmt 
Thus, re aS oan) Fe m*. e 


The equation may now be written : 
m?, emt 3. 7m emt + 10 emt = 0 
Dividing by e™* [emt + 0] 
m+ 7m + 10 = 0 
“im = 25)or =2 


y Xe-5t + Yer-2t 
To Check. 
iy = =/5 Ke — 2 Ye°* 
ek 4 
al * = + 25 Xe-st-+ 4 Ye-% 
LI ed st a 
= = ats ~5t 
2 a er + 10y.= 25 Xe-*t + 4 Ye 35 Xe 
— 14 Ye-** + 10 Xe-* + 10 Ye-2 
= 0 
Consider an equation of the second order, 
ay SAAT ear : 
EEN Tea 
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This may be written : 


m+n? = 0 
m= — 0 4 
m= tin [FJ =vV-T) 


Aont+ Bevint 

= A[cos nt + 7 sin nt] + B [cos nt —7 sin nt] 
= (A+B) cos nt + 7 (A-B) sin nt 

C cos nt + D sin nt 


Il 


y 


Application to a Mechanical Pendulum. 


From Newton’s Second Law of Motion, F = d (mv) / dt 


dv i dm 
= = its 
m di u a 1 


If the mass is constant, the second member of the right-hand 
equation is zero. 
Force = mass x acceleration. 


It is usual in examples of this kind to use a “‘dot’’ notation for 
differentiation. 3 
: if « = displacement 


Fig. 2.—Mechanical Pendulum (Newton’s second law)} 
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velocity v= 


acceleration = 


Referring to Fig. 2: 


T cos 6 = 
- Tsin 6 = 
length of arc = 


Thus the Restoring Force 
= /sin 0 

& 

tan @ 


dx = 
57 = @ 


dt 
ax 
dt? 


mE 
m & 
16 


= mlé 
1 @ [@ is assumed small] 
= 16 : 


(-1g) 


! or et tan 0 


6 
0 
6+ glo =0 


This may be compared with : 
d*y 


The solution is, with the premise, 
6 = a, 


6 =0 when ¢=0 


@=acost V gil 


Alternative method, 
6 is not assumed small. 


Fig. 3 is an elaboration of Fig. 2. 


It is known that potential energy + kinetic energy = constant. 


involving energy considerations, where 


P.E. = gml (1-cos 8) 


K.E. = $m (il 6)? 
gml (1-cos 8) + tmP 


Differentiate to eliminate K. 


gmlsin 0+ 4mP 2 


= gmilsin 0 + mP6= 0 
6 + g/l sin 0 = (0) 


@= K 
dO d ao 
(ay ae PH pean 


(as previous solution). 


mgl (1—cos 6) + 4m (10)? =K 


(8)? ='2n2 cos 0 + K [n 


= Veil 


Fig. 3.—Mechanical Pendulum (Energy method) . 


Let a be angular extent of swing, 


=) 
> 
=] 
a 
= 
i] 


[= 8] 


a when 6 =0 


= 2? cosa + K 


-2n® cos a 


2n* [cos @ — cos a] 


= 1 - 2 sin? (a/2) 


1 - 2 sin® (6/2) 


fu a0 
J 2 4/[sin? (a/2) — sin® (6/2) ] 


1 
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The solution of this involves the use of elliptic integrals* and 
the time of a coniplete oscillation is given by 


1 a 
T=4.—.F]si > 
I 4 7 sin > | 


F [sin a/2] is the elliptic integral of the first kind (modulus 
sin a/2). 
Expressing as an infinite series : 


27 120s a 8 d a 
ar a [1 +o sin® (5) iy Bau ye Se (5) 


120032652 ren /Acl 
+a pen ()+---| 
A first approximation gives : 
T= 2itin = 2 V Ile 


Comparisons between Electrical and Mechanical Systems. 


ro 


ad a 


Resistance. 


(a) ELEcTRICAL. 


Effect: to cause dissipation. 
R = eft 
where R = resistance (abohms). 
e = voltage across resistance (abvolts). 
4 = current flowing through circuit (abamps). 


(b) MeEcuanicat (Friction). 


Effect : to cause-dissipation of energy. 


P = flv 

p = mechanical resistance. 
f = applied force (dynes). 
v= 


velocity (cms./sec.). 


* SAlaV da] da A is a rational function. 
¢ (a) is an algebraic function. 
Integral of the first kind. 


fae da 
F (8, «) = fe (=e) (I=a? a) 
ay 
S=8 d8 
= | V/1—o sin? & 
an 


Provided the assumption a = sin-} 8, and w? < 1, holds. 
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Inductance. 
(a) ELECTRICAL. 
Effect : to oppose any change in current. 
on 
dt 


voltage across inductance (abvolts). 
Inductance (abhenries). 


eé= 


iol 


é 
Jey 
di i 
ar Rate of change of 7 (abamps./sec.). 


(b) MrcHANIcAL (Inertia). 
Effect : to oppose a change in velocity. 
dv 


af PLO! oy 


m mass (grams). 


iol 


fi applied force (dynes). 
& = rate of change of velocity (cms./sec./sec.). 


Capacitance. 
(a) ELECTRICAL. 
Effect : to oppose a change in voltage. 


° . dé 
t=C. ‘a 

C = capacitance (abfarads). 

7 = current (abamps.). 

. = rate of change of voltage (abvolts./sec.). 

F : 
‘le 
e= =| 4.dt = 5 
where g (=fi.dé) = charge on condenser... .. abscoulombs. 


(b) Mrcwanicat. (Stiffness Compliance). 


f-Z-s 


compliance (cms./dyne). 
1/C. stiffness (dynes/cms.) 
applied force - - (dynes). 
displacement (cms.). 


RSNO 
hoki 


' 
1 
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DEGREES OF FREEDOM IN SYSTEMS. 


In Fig. 4 a mechanical and electrical ES yetem is shown, each 
possessing one degree of freedom. 


C 


R 


Fig. 4.—To illustrate single degree of freedom in mechanical and 
electrical systems. 


Fig. 5 shows two similar systems, each possessing two degrees 
of freedom. 


- ! 5 52 
| | Pe 
| 
; | 
| 
~ 
t 
Me 


— Xr 


Fig. 5.—Double degrees of freedom. 
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Equations to be used in the Analysis of Systems. 
MECHANICAL : Newton’s Laws. 


e.g., m% = applied force — Restoring and frictional forces. 
I @ = applied torque—Restoring torque and Frictional 
forces. 


ELEctTRICAL: Kirchoff’s Laws. 
e.g., applied forces (e.m.f.’s) = 2 back (e.m.f.’s). 


PARTS OF A DIFFERENTIAL EQUATION. 


The equation, am?+ bm +c = 0 is the Auxiliary Equation. 


A solution of the form, y = Ae**+ Be-™*, where A and B are 
arbitary constants, is the Complementary Function. é 


When values are substituted for these arbitary constants, the 
solution is the Particular Integral. 


BREAKDOWN OF EQUATIONS. 


Consider the equation, 
; by 
sr +6 +9y = 0 
dt dee ie 
This reduces to: m?> + 6m +9=0 
(m + 3)? = 0 
y = Ae-%*+ Be-% 
= [A+B] e-*= Ke-* 

Wherever the roots of the equation are equal, the equation must 

be reduced to the form, y=e™ [Ax +B]. : 


Consider the equation : 


dy dy Da eepe 
7. + a.a t b.y = Ke 


If B is a root of the auxiliary equation used in obtaining the 
complementary function, then the particular integral cannot be 


found by making the substitution y=A eft. 
The substitution, y =At eP* must be made. 
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APPLICATION OF THE CALCULUS TO BENDING PROBLEMS. 


The following relationship is derived from the theory of simple 
bending : 


1 MM ss ole a 


Ig) jaa ii y 
1 dy : ue 
eae de (when dy/dx is negligible). 
dy M ay 
TI = P = WM 
hus, ea ai or EI de M 
d 1 
Also, = = EL .w (Loading) 
dy 1 r 4 ‘ A 
= soe | ow. dx (Shearing force (S.F.)] 
ay 1 fir r 
a = Ge | | w . dx . dx [Bending moment (M)] 
dy 1 P 
rele ee =r | | | w.dx .dx . dx (Gradient]. 
1 , 
G0) = a | [| | o-er. ar. ar. dx 
EI JJJ) (Deflection] 


(1) With reference to Fig. 6, 


M = = (L-x)2 
ay w Se 
PR op Merde 
dy w 
Bt en eis 
in 6 (L-x)? + A 
When x =2- ee 0 
ax 
2 3 
(ye i ee) ee 
6 6 
dy w wl 
os ==) 
Be Ba an tel we 


AN 
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M. a 


RO" 


Fig. 6.—Cantilever Beam, uniform load. 


w : wBx 
Ae oy Uncen ee 
when x=0, -v=0. 
wl w lt 
ns B= 
Os pated 24 
: w w Bx wl 
aia = Seay pr 
Ely aa | x)t + 6 4 
1 wilt w 14 
EI y (at x=l) aT wl 
= wl ay 
— 8 aS 
wl = W 
Wi 
2 SE 
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(2) With reference to Fig. 7, 


d*y w nel 1 
M = EI we = oy (ae)) - a C-*). 
Fe ee ea A ai a A : 
dx 
lw a 
Bie = 2 G4) Seep 
2 m4 ¢ x) 1 | ae + Ax +B 
when x = 0, y=0 
x =1, y= . 
wl4 w 14 , 
EI y (at = = an de 13 28 0) 
y (at x = 0) wl Tone 


Fig. 7.—Simple Beam, uniform load. 
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(3) ‘Reference to Fig. 8.- 
Fe ; M (at ®) = De 
ay W 
EI = 
dx yaa 
4 Wie 0% 
y (at x = 1/2) = = 
ee ate 
, um 
| jetee i 
Tanne) 
, eee 
ip bea |e ae Ss LEE 
rj ‘ 


. = 
‘Fig. 8.—Simple Beam, concentrated load at centre. 


PARTIAL DIFFERENTIAL EQUATIONS. : | 
an problems pertaining to wave-motion, the following equations | 
exist :-— : : 

y = f (%—vt) representing motion to the right. 
y = F (x + vt) representing motion to the left. d 
The equation y = f (x—vt) + F (v+vt) represents the double 
wave motion. i 


22 Ff (eo) + FY (e+ ot 


oy 
= —uf («—vl) + uF’ (x + vt) 
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ay ‘ 
: = f” (x—vt) + F" (x + vf 

ae f” (% —vt) (x + vt) - 
G 

os = v? f” (w—vt) + v? F” (x + vt) 


(1) 


ay 1 pay | 
Ox? v2. oe 


Application of above Equation. s 


Considering a section of an ‘“‘ideal’”’ transmission line having a 
capacitance of C. farads per unit length, and an inductance L. 
henries per unit length. On a section of the line, 5x, in length, 


Current = C Ow 
Gt Ox. 
de ot 
Moltacms= namie ao 
oltage oe ry 
Differentiating, we have: 
3 CBE @e 
ie = a {with respect to ¢] 
Pe ib Gee ith ic ite x] 
= respect to 4 
ax Pre vines : 
@e Be 
Thus, = IL 2) 
0x? LE 6 


By comparing equations (1) and (2) we have : 
1 
a Sole G; Wes 
or VLC 
The capacitance of two parallel feeders (air dielectric assumed) 


Salle A 
BSS See EG . 
Aeh aly arads/cm 


(where d = spacing in cms. and y = radius in cms.) 


iw 
a 
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Employing the same symbols : 
Inductance, L = logh d/r x 10-® henries/cm. 
1 ; 


UC 


Vi = 


3 ‘ fee . ——— . 10°.9.10-2 2 
; logh d/r Tea 1035-29) Jo 


phi a 
Y 9 x 102 


= 9x 10% = 3 x 10% cms./sec. 
This is, of course, approximately equal to the velocity of light. 


Thus waves travelling along air-surrqunded lines attain a velocity 
approaching 186,000 miles/sec. 


INTRODUCTION TO HEAVISIDE’S OPERATIONAL 
CALCULUS. 


Contrary i common opinion, Heaviside’s operational calculus 
is not a new branch of mathematics. | Similar methods were used 
by earlier mathematicians ; its originality lies in the application. 

It has already been seen that the first derivative, or differential 
coefficient, may be of the form, dy/dt: 

This could be written d/dt . y. 

* Heaviside writes for this; p.y; thus p is equivalent to aie 
and it is an operator of differentiation. 


Similarly the symbol 1 /p is the reverse of the symbol p. There-, 
fore 1/p is the reverse of differentiation, namely integration. 


d 
Thus, P= a 
1 is 
ie ie ].dt 
Example. ? 
To differentiate 3x2 + 2x +14 
i.e., p (3x? + 2x8 + x4] = [6x + 6x2 + 4x3] 


orp. 3x2 + p. 2x3 + pint = 6% + Gx? + 4x8 
Thus, the operator p may be used algebraically. 
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Example. 
Differentiate x5. 
be) = put = 5x! 
Differentiating again, p (5x) 
Differentiating again, p (20x3) = 60x, etc. 


i] 
to 
oO 
8 

eo 


n = 


dit 


Thus, p? = 


; Example. 
Integrate 1, with respect to 4. 


‘Thus, ai =-t ‘(- +) ; 
BAAS 1 

nteera cine aa Liesl £ 

ntegrating again, = — = = — 

grating ag: 2, 1.2 [2 

Int ti : tT i? te ie 

ntegrating again Ss lire | 123 743 
*: 1 B t4 14 

Integrating again = + [ = == 

” PP Lb23. 1234 4 
Thus, alee ua 


B ~ [a - 


A cautionary note must be stressed, As the symbol / is an 
operator, it only operates on terms following its insertion.in an 


‘equation. Thus x? is an operation demanding the differentiation 
of.x?, but x* is not. ; 


THE UNIT FUNGTION. 


Suppose that a force is about to be applied suddenly to a system; 
before time ¢=0, the magnitude of the force is zero, but immediately 
subsequent to ¢=0, the force rises instantaneously to a specific 
value, and remains at that value. Thus with such a supposition, 
many complicated problems may be reduced to comparatively 
simple calculations. Heaviside visualised such a condition, and 
based his operational calculus on this. The specific value of this 
force or impulse, being unity and represented by the symbol (1) 
See Fig. 9. ' J y 


FY 


~<— (t-}) —t+—(t+fa) — a 


t=O 
Fig. 9.—Representation of Unit function. 
{1] =Oatt <0. . 
(lf =latt> 0. ; We ac 
‘ ' [1] = (toljatt=—0. o 


An example will now be demonstrated in full, by three different 


methods. 


Example.—A “trip coil” is energised from an uni-directional 


current source. See Fig. 10. 


s Tere (Gm 


Fig. 10.—Elementary circuit of trip-cdil energised from a battery. 


(1) _By THE ORDINARY CALCULUS. 


It is known that : ‘ J Z 
y Ri + Ldi/dt = E ~ (3) 


Adopting the usual procedure for eae differential equations. 
Ri + Ldijdt = 


Let ‘“D = djdt 
then, (R - +LD)i=0 
Let 7 be of the form A e™ : 
“then, difdt = Ame™t : 
a AR.«™ + AL.m.e™ =0 


but e™t + O 
- AR+ ALm =0 ¥ 
r . ALm = - AR 
: m = — RIL 


iy = Aé-RIL-t. 


, 
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The steady-state solution of the circuit is obviously is = E/R 
Therefore, the full solution is, 7; + is. 

‘e 


Hay Ae-Rh.t. 


If the time is measured from the instant the circuit is completed, 
then 7=Q when ¢=0., 


Ig] deh Ware = 4 Ae-Riket. 
: 0 -= + Ae? 
A jes we. 
| R 
| = E oe E e-RIL.t. ne 
m7 


ile = | 
\ 
(2) By HEAvistpE’s METHOD, LEADING TO AN INFINITE SERIES. 
! Ri+ Ldijdt = E h 
| OR. ab A) eS a (4) 
£ te en (4a) 
R+Lp 
E ‘ 
ee eee 1] (5) 
LP L(R/Lp) +1]! 


ral 


let a = R/L 
- oe 
— Lp l(a/p) + 1 
1 ; 1 
ip?) Pp apy Ol a 


SI : 
Evaluating [ + al by the Binomial Theorem. 


at 


} = 2 2 
a ey feet GP ok tip 
Se ies peat wt pe 


i== 5 es a 
at a® i? ak & at t4 
--Sb-n* Bp ere pete. |] 01 
Be goelp, ee Be qa 
L [2 [3 14 [5 (7) 
Ifa = R/L, then La = R and L = R/a 
J ta cee atts a3 tt atts ' 
then 7 = zr ee + [3 = [4 + IE te | 
E att at ttt atee 
ly Let Tetsu. ae Ma tc. |, 
Referring to the Infinite Series for e +x jt is known that : 
% x x x4 x 
ie mama eee) (1G) 4" ~ |n 
at BRI $2ERE 20S aa 
oF Sea mias 8) le ac 
Therefore, substituting at for ‘x, it is seen that : 
\ at at? a’ 8 atts 
gout 29S a at we - [3 + 1 etc: (10) 
at at? aS fs a4 74 
+ + etc. (11) 


Se Rady coma Sher 2 has 


To obtain equation (11) in a form comparable with equation (9), 
multiply throughout by (-1). 


re Seer es [i [2 [3 = [4 
4 = B/R [Pl — es] 
But a = R/L 
BS IDYIR, [Lil estbeots || 


i] 


The above method may seem somewhat involved, but a simpli- 
fication may be made by using the Expansion Theorem. 
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THE EXPANSION THEOREM. 


E 
In equation (4a) the expression ¢ = ——— [1 
q (4a) xpression 7 Raap |! ] 


was deduced. 


z.e!, The current i = The voltage E, divided by an “impedance 
operator,” R+L%, operated on by unit 
function [1]. 


The operational solution of this may be written : 
j= E22 (12) 
So) 


ae $i) and fip) are functions of p, fp) may be a polynomial 
in p, viz.: 


Se) = p + apr-2+ B pre + apr-e + Spn-4 xpr-* (13) 


where x is the (w + 1)th term (obviously a, }, c, d, etc., form an A.P. 
series whose common difference is 1). Figures are not included, 
so that the series will be as general as possible). 


The roots corresponding to f(,) = 0 are assumed to be : 
(ay the HG - - - ita 
“fey = (P=) (6-72) (P -75) (p —T) (14) 
dp) May also be a polynomial in #. 


«ffi is of a higher power in p than ¢), then obviously ipy/fip) 
will be fractional. 


Conversely, if ¢() is of a higher power than fi), then by division 
Pi/fip) Will result in-a series of terms containing one fraction 
concerning the two functions. 


; From the foregoing, it is possible to write : 


Pe) ps So) 


Poy ah SEP ADL ew 1 AL Oh 8 
fio) (P=%) (P=72) (b=7) etc. (15) 
Solving by Partial Fractions : 
$¢0) A B Cc 
WAS (se 16 
Fa (PAS. Wl. he 
\: to. oo. . A f Do 7) 


Foy a (p-1) Fo (2-1) Fo 
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Multiplying by (p-17), 


dw) Poy 

Se eS aay) 18 

Fo mob Ce 
when p = 7% 

Pp=n) 

Sad aN 19 

Fe=n) ( ) 


From equation (17), 
fe) = (-n) Foe {comparable with a product : 


(aera dv du 
See) aaa ae 2 || 
fie = Fe + (6-1) Fw (20) 
when p = 7‘ 
fy) = Fo 
e Mestaly Vie (21) 
fas 
Similarly : 
Se) =B (21 a) 
fe) 
Se) =C (21 b) 
fe) / 


and similarly for each root. 
It will be appreciated that: 


fe) = B - fo) (22) 


and that fu Fa f(y ete., are the corresponding values of f"(p) 
when the respective roots are substituted. 


Thus equation (12) may now be written : 
5 A B Cc 
t+=E [ + + 
| P-n p-'s P-Ts 
It has already been proved that : 


ier alse 


+ etc. |a (23) 


weal es 5 
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Substituting —a for +a 


pra tp] 


This may be compared with equation (23). 


peB[ -2-3 -f-. ete Jay 
% To 13 : 
Ae B erst C etst 
+ [ + + ere eLC., ] [1] 
n Te T3 ‘ (24) 
The first paranthesis may be compared with : 
A B & for the condi- 
[ + + He tehsteyie | tion p=0. 
p-7, P-1re P-Trs 
[-4-8 -f-...te ]-[] os 
fi v2 13 fey 4 p=0 
which may be written bo 
fo 
A= [ fo] when p=7, may be written A = Peay 
fo Ses) 
In the same way, 
B= Pees etc. 
Ses 
idl pe $0) | if the condition #=7, obtains (26) 
% Pf’ 2 
Thus it is possible to write : 
bio) Pt lean | 
+etc.; [1 27 
# dflap pafiep J, * ey 0 ee 
Thus equation (24) becomes : 
%o , yoo al 
115 le (1) 28 
to ra (P) es) 
Tey UE Ue f= 


This is known as the expansion formula. 
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There are various methods of writing this, one of the most 
common being : 


eee (ho) pow 
fe [ i a | a) (29) 
Zoo) *P aZldp 
pp ji Basan 
Application of the Expansion Theorem to the Example under 
Discussion. 
1 
pS by ell 30 
t= ER ap ag (30) 


Re-writing in terms of equation (29), and arranging in a form 
conducive to application of expansion theorem, 
Cera se (LU (31) 


Thus Yoo =1: Zp =R+1p 
Equating Zip) to zero, one root is obtained, viz. : 


R 
ih SS it, 
dz dz 
= Land p— = 
dp ne? dp i 
eo = pa (32) 
p dzldp pL 
when p=); 
845) ly SEZ OO salt (33) 
p dzfdp I_, -R/L.L R 
1 
and Yo (i.e, p=0) = R 
Zo) 
G Yi, pt 
Pe [s+ 2 pS aemiap na (34) 
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ROOTS. 

In the development of the Expansion Theorem, the roots 
1%, Yo, 3. » - Yn, Were extracted when f(p) (this may now be termed 
Zw) =9 The treatment accorded to these roots obviously 
assumed that: 7, 2 7% 273... 2 %, but that: 7 + 7% + 
13... +%, +0 and that the powers of 7%, 79, 72... Yn, are 


integers. Depending upon the type of equation formed, many 
interesting variations of roots may be investigated. 


EVALUATION OF A POLYNOMIAL 


Due to Horner (1819). 


Supposing a value of f (s) be required for specific values of s. 
The method described below will be found to be extremely compre- 
hensive. j 

Assume a polynomial in s of degree n=3 : 

asé+ Bs? + 8s+60=0 

Arrange the coefficients as follows :— 


ny we 3 0 Ls 
as as*+ Bs ast+ Bs*+ 8s 
as+B as*+Bs+68 as? +Ps?+8s+0 = f (s) 


The operation will be obvious; namely, multiply the first 
coefficient (i.e., the coefficient of the highest degree of the function), 
which is placed as a leading entry, by the particular value whose 
function it is required to find. The product as being placed 
underneath the next coefficient, and added. The result is again 
multiplied by s, and added to the following coefficient and so on. 

A method of synthetic division may be found extremely useful 
for the division of, say, f (y) by (y-a). 

Let f (y) be of the form: ay®+by?+cy+d = 0 


PROCEDURE : 
a +b +¢ +d (a) 
+aa + (a’a + xb) + (a8a + a®b + ac) 


a@ + (aa+b) + (a2a+ab+c) +(a8a+a*b + ac + d) 
The quotient is: ay?+ (aa + b)y + (aa +ab +c) 
with a remainder- a’a + a*% + ac + d. 
Determination of Irrational Roots (Horner’s Method). 
f (a) = Aat + Ba? + Ca?+ Da+ E =0 
Reduce the equation by the elimination of rational roots, thus 
assume (see preceding method of synthetic division) the equation 


becomes : 
Ba@+Da+E=0 
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By various methods it is possible to find the limits affecting 
any change in sign in f (a), thus it is possible to find a root which 
lies between 0 and a. 


Dividing f (a) by (a-a): 


B +0 +D +E (a 
+aB + a®B +aB+aD 

B +aB +@B+D |+@B+aD+E 
+aB + 2a°B 

B +2aB +3a7B+D 
+ aB | 

B +3aB 


The transformed equation now becomes : 


Ba’ + (3a B) a?+ (3a2B + D)a + (a3 B + aD + E) =0, 
and possesses a root between 0 anda. A more exact approxima- 
tion of the root may now be found. By diminishing the irrational 
root of F (2) = 0 by, say, an amount 8 (where 6 is smaller than a), 
a second transformed equation may be found. 


B +3aB_ + (3a°B + D) + (a®B +aD +E) (8 
+ 6B +8B+368a + 8[(3a°B+D) +(8B +38a)] 


B +B (3+3a) + [(3a°B +D) +(5*B +38a)] | +8 ((302B +D) 
|+ (8B +38a) ] 
|+ (2B +aD +E). 

+ 8B + 8(25B +3aB) | 


B +(28B+3aB) |+8 (28B +3aB) +[(3a2B +D) +(52B +38a)} 
+ 5B | 
B+ 38B+3aB 


This equation resulting from the second transformation will 
be found to possess a root between much narrower limits than 
previously. Thus the root of the original equation is : 


a+S8+A+7+¢ + etc, 


where A, #, ¢. . . etc., are roots of subsequent successive trans- 
formations. 
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EVALUATION OF A POLYNOMIAL BY DETERMINANTS. 


Assume the following simultaneous equation : 


a,a+ Bb = 5 (36) 
a, a + Bod = 55 (37) 
where a and 6 are unknown. 
Thus: . a,B,¢°+ B, Bob = Box (36a) 
a2 Ppa + B,B2b = By x» (37a) 
(a By @ — a2 ByQ = (By %— By %2) 
cis ta fee Box — BX» 
ay By — ay By 
’ Adopting the method of Determinants : 
ie - By | 
ee (altar Bo (% By - %2 B) 
ay. By (ay Bo oe B:) 
ay - By 


The method is quite obvious. N; amely, coefficients downwards 
from left to right, constitute the positive terms (depending of course 
upon their individual signs) and coefficients upwards from left to 
right constitute the negative terms. 


GENERAL CASE FOR THIRD ORDER DETERMINANTS. 
| a. By. & a, . By 
AK =| a. By. 8 a, - By 


as - Bs - 83 |. a; . Bz 


The terms a, f, etc., are written to the tight of the square array 
of elements, merely for convenience. 


Thus : 
ae kc tae a Bo dy 
> Bs 82 a4 — 83 ae Pi] 
[oy (Bo 33 — Bs 83) — G2 (Bi 83 — B; 3,) ; 
+ as (B; 8 — Bz &) ] 
ay By ey) By 8. | Br 8} | Br & 
a, Pp 8» = a | + a3 | Bo 8, 


| | | 
[as Bs 8s| | | Bo 8s 


or, 
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The lower-order determinants are found by suppressing the 
column and row corresponding to the particular minor. 
A = a [Bz 83 — Bs 85] — a2 [B, 83— Bs 8) + ay 
[Bs 8. - Bs 31] 


TYPES OF ROOTS, BY CHANGES IN SIGN. 


It can be shown (e.g., Lill’s graphical method) that there will 
not be more positive roots than changes in sign of f (a) ; similarly 
for negative roots with f (-a). 


Example. ; 3 
f (a) = a+ a8-5a?+ 8A-9a+1=0 


In this equation there are four changes of sign, so four positive 

roots may be expected. 
f (-a) = a+ a8-5a7+ 8a-9a+1 =0 

In this, there are two changes of sign, corresponding to two 
negative roots. ~ 

Actually, of course, in the example discussed, there will be ten 
roots. Six of these have been accounted for, and it is to be ex- 
pected that the other four will be imaginary. 


Let f (k) be of the form : 
a (ka) (kB) (k-8) 
It may be that a= or B=6 or a=5, etc. 
Obviously (k — a) is a factor of f (k) if a is a A-repeated root. 
f(b) = @ (R-a)A (R-B)o (&-8)4 


If the equation is a n’th degree equation, then A+ ¢ + + 
succeeding powers = 2; and all derivatives up to and including 
the one preceding the \ term will cease to have any significance 
when k=). 


( 
( (e+ $ (B) A (b=a)a-? 
f" (&) = (k-a)\ 6” (k) + $' (A) A (R—a)A-1+ g (k) 
A (A=1) (k-a)A-®+ ’ (Rk) A (R-a)A-? 

It can be seen that the first derivative contains the factor 
(k—a), (A—1) times; the second derivative (A—2) times ; and it is 
clear that successive differentiation would give (A—3) times for 
the third derivative, until finally the \’th derivative would eliminate 
the factor. This is the procedure to be adopted when repeated, 
or equal roots are present. 5 : 

When one or more roots are equal to zero, or some roots are 
repeated, it is not possible to use the expansion theorem in exactly 
the same manner as previously outlined. : 


. 


FE 
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It is necessary in examples of this nature to perform an inte- 
gration. 

Thus,if *#=E.Yp/Zp [1] 
where Z) is of the form, 


P (P-n) (b-72) (b-713), 


then obviously one root is zero. 
The procedure would then be : 


: Yop) 
Evaluate 1 = E . ——— LSESEEEN ES (1) 

i (b-%) (6-72) (6-70) 1p 
then perform the operation, which is demanded by 1/p, (the missing 
term in the equation for zo). 

The integration to dt must be performed. 
a0, 
When more than one root is equal, a “‘dodge’’ may be resorted 
to. Each root is imagined to be slightly different. 
Thus: [p — (4, -8)] * [6 - (4 +8) J 
As these two approach equality, the necessary corrections are 
made. 


QUADRATIC EQUATIONS. 
_ The equation, ax?+bx +c =0, is a quadratic equation, 
1.€., a polynomial in x, of the second degree. 
The procedure is as follows: 


Add [b/2a]? to each side. 


2 2 
e+ 2x4 [2]- -£4[F] 
a 2a a 2a 


P[esf] Sete 
2a. 4a? 
E i b 2 ey B = 4ac 
2a 2a 
-bivV B= 4 ac 
WT ae 
; where a +0 
For equal roots, 62 = 4ac. 
For real roots, b> 4ac. 
For imaginary roots, b? <4 ac. 
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The quantity (62-4 ac) is known as the discriminant of the 
equation. 


Let b?- 4 ac = wo 


If wp) > 0, roots are real and different. 

If wp = 0, roots are real and equal. 

If w < 0, roots are complex and different. 
If wp is a perfect square, s?> 0 


-(6+4s 
then + = eel EE and both solutions are rational. 


Whenever complex roots of a real equation occur, they occur 
in pairs, e.g., a +ib. These are conjugate complex numbers. 
Obviously the number of complex roots must always be even. 


The following problem will be discussed at length, as several 
important principles are involved. 


Problem.—A resistance condenser and inductance are connected 
in series across a steady e.m.f., e. 
The voltage equation is :— 


ey wey ee OL | aa =D (38) 


Re-arranging and differentiating : 


: Li : 
Ly pA egy en al (39) 
1G dt dt TAG. 

Adopting the symbol D = s 
then : 
R 1 : 
p+ E.p+ —]i-0 (40) 
[ It; ic 
- R/L + V RIL? - 
raed | as | 4/LC (41) 
] 2 
(ial oe 
Dp DEEN AT eat ya Eas (41a) 
2L 41? IEG 


2L V ae LC 2L Va Lc 
eee + Be (42) 


. (-=+ peReAna la (-=- [Say 
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+ - aba (42a) 
=a |g = 2a 
Wes [ my) ae area 
x Vy pak (42b) 
tlle ice ier ELC 
i = Aet + Ber (43) 
ane t = : be of three f 
e term Lae ‘ : 
412 LC may be of three forms 
R? 1 
(1) > — 
4L? Le 
R? 1 
(2) a 
41? nG 
@) ae 
4L* LG 
Goes (iy Na 
Zt ILC. 


From equation (38), if at the instant of switching-in, the con- 
denser is not charged and the current is zero, then : 


di A 
— = 1D. 
dt 

di E 
eae Ts 


Thus equation (42) becomes, 
Ae + Bel = 0 = AFB 


It is possible to write, 


di 
9 = pA + p, Bert 
If ¢=0 
dt 
= pA + B 
a Ar Pe 
di E 
But —— = — 
ie ve L 
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pA + p2B = E/L 
Bo eae - 


PA - pA = B/L 
A (fi- #2) = EL 


Q 1 
Ac aes 
L (pi ~ bs) 
1 
endo e 
L (p:-p2) 
Equation (43) may be now written : 
pe B 1h RSE ae 
L (pipe L (b:— Pe 
E 1 nen 
ath [emt — ert] 
L (Pb: - ds) 
R Ta 
B = Poet i ae 
Me aehe OL 4L2 LC 
R R? 1 
De) SO en arate 
20 4L L¢ 
E e~Ri/2b 
i a 
2 
Tay) es ee 
41? LC 
my =, 
[. (Jam -ae GA an a) ] 
—- é 
eb — rd 
It is known that, sinh 6 =" 9 
2 
P 1D es eel iy Re ee 
Thus, 7 ( en —)+ | 
= es sa ERECT 2 Tm gol Sinh ren Le 
412 


This is an overdamped condition. 
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Case (1)—By HEAVISIDE’s OPERATIONAL CALCULUS. 
Equation (38) is now written, 
ee a1. hi 


ae pic = fii 


Thus : 


ro) 
Hoon 


To obtain roots of Zp) = 0 
Zp = LCf#? + RCP+1 =0 


a R / R 1 

2L Ay Th we 

Thus : 

R R? 1 
bh = -—+V— - 

oe Ay esTG 
4 a iP ve 
i 2L Ane 7 Te 


Zw = (2LCH + RC», 
= tvVRC? = 4LC 


a) Saat uel zak 


41? LC 


The expansion theorem from equation (29), 


xf z[ 2 Yo Yo 7) iy 
ay  y|dp | 


Zo) Pu »Pay Pa 

0 Cpe 
we =E oe + 2 R 7 ] 
} pom p 210 ne 


Be ta Le nT ioe” M x A Sk a i dhl a ae ol 
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Pl ert 2 pat 
- Ec| 0 5: =a ee 
2 
pr 2LC ee pie: yeas 
412. LC 7s ie 
Rin 
Let B = (= - =a 
412 LC 


*Y 
= EC[o + Brae aly as ere 
Pi2LCB p22L C B. 


¥ E.G py ePat po ePst ; 
acpi san heal 


ye [emt — erat] 
2LB 


E . en RAL 


rape (ree het 
412 = Le é ic 


From the previous definition of sinh. $, the above equation 
may be reduced to: 


2L 


-R/2L. —!—— += 
Ngee dee A R? 1 
iy) ee [ sins ( ee)" || 
5 Nee ie G 
Case (2). Re 1 
pep nk 
412 iD 


In this case, 1/LC is greater than R*/4 L*, thus the quantity 
inside the square root sign is negative, therefore the roots p, and po 
are imaginary. 


f 1 R? 

ue abicnin iene Se 

) A Apis cme a 
5 20 AE TZ 


(In this particular example the symbol j = V1 is used in 
place of i = V —1, as ambiguity may result, by reason of the 
current symbol 2). 
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‘Sulbsiaianiaing in the coe aa 
7 A ePit + 3B est 


we have 
fie Reaves fa RP R and =] 
G 2L Vere earl” [- aL 7VLe ate 
aa Ale a + Be 
Be ijlme RURe - eee 
[ Times aalk ne os ] 
= e RPL+L Ae + Be 
8 + id 
cos = ——— 
2 
peed x j8 — ¢-i8 
similarly, sin @ = eae 
23 
20 + eid [eid — eid ei8 + ¢i0 20 — 6-16 
Daa al ] Rhee peed 
= 6p 
<% 26 = cos 0 + 7 sin 6 
also, 
e38 = cos 0 — j sin 6 


Therefore, the expression reduces to : 


Gls emattl A( cos — = x +7 sin me Le =) 


4 
R? 
(VES jenVE 2) 
(cos ane —j sin ane 
1 R? 
reser “Pasay ) cos eee _~— j (A —B) sin 
Te > funy elk ) 


1 4 Re | 
1G es Oe 


URE Be / 2 
then, t=e8?L-# 4/C2_D? sin ( Vv = = .t+¢ ) 


Let A +B =C and j (A-B) =D, 


Let M = VC =D? 


ie 
3 = is 
i= Me-R? + sin (ye ~ Ft +4) 


oO 
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Fig. a URS to R.L.C. circuit (case 2). 
With reference to Fig. 11, 
tan ¢ = C/D 
or, ¢ = arc tan C/D. 
With conditions obtaining as in the previous case, 
Me°’sin 6 =0 (M +0) 


R? 
Thus 7 = Me-R2L-* sin an t 
LC Ay? 
«ei ay 
ee if oan ee Ne. Zune 


If ¢=0, then sin 0=0, and cos 0=1. 


di i V. TER? 
- [m. -=,.1-0 | 
di LC 41? 


4 di 1B: 
Buty =. == 
an a L- (see Case (1) 
2 
Re ie ae 
LC 4L? ib, 
M = E. 
4 A. 
LC 41? | 


ws ot eh - md OE ne 
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The full expression now becomes : 
2 1 R2 
aa i DS Ceo (sin vt -_— 
Lge Le Pia! 
LG 4° 
This is an underdamped condition. 
Gase (2). By HEAvIsIDE’s OPERATIONAL CALCULUS. 
Ri+Lpi + yP.ifc.i=E [1] 
al 
a Rue [1] 
R +L + 1/P.1/C. 
Thus Yq) = 1 
Z eRe Aes 1P~ 1/G. } 
i Bey, = 0, 
then, ith 
py - R/2L +7 V1/LC — R741? 
pb, = -— R/2L - 7 VIL — REALE 
Zp) = L — 1/p?C 
- PZ") = pL - 1/pC 
LF RILG 
pc 
2 
} [ez | 5) py See 
pi IE (G; 4L? 
Tay aa eee 
[2 Z' | = -2Lj Pah IN 
Po Ie) Ya 
Yo) = 0 
Yoo) fe 
Zo) 


Substituting in the expansion formula : 


Past Vecc acl Pawel 


Poe 19; (0) tote ie feet 
2 
2Lj (ae ek a * 
1 Le 4L2 a 
& [1 [fay 
(vee aah ohn ie a) 
E e-R/2L.t eee 
*eeeepal 2 
L 1 R J 
1G 4L? 


wee ee et Oe vo re mes 
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Sin @ has already been defined as : 


cid — id 
27 2 
F Premises : ma i Ri 
ye LC zt) ] 
ie 4L? 
| Case (3). co 
Mes Tie 
R? 1 : R? 1 
: I = —— the expression —_ - — 
4L? LC 4L* LC 
is equal to zero, therefore the roots are equal to —R/2L. 
7 = A e7R2L.t = Be-R2L-# 
= (A + B) e-R2L-t 
1 . 
, Thus (p* + Dy, + —) zt =0 
It, LC 
2 2 1 
[ (P + R ) Re 7=0 
20 412 LC 
5 f R? i 
t = —— , 
Bat Ti Wa ae 
pines ji=o 
( v ORey 


\ R : 
Let A = [p a | 

ue 

then 0 = 4[D + ea 


From preceding discussions, A = (NGS 
Thus) | (DeeiR/2 ye — Aten Sat 
(D + R/2L) 7 Aapiaol IN 


as before, D = d/dt 
, 4 eRRL-t = At + C 


where C is an integration constant. 
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Multiplying by e-®/?4-+ 
3 SUNG 4 OC) Geese 
The boundary conditions are the same as in the preceding cases. 
E/L=A and C=0 
Thus the final expression is 
4 = B/Lot.esR2t-t 


This is a critically damped condition. 


Case (3). By HEAvVISIDE’s OPERATIONAL CALCULUs. 


Rearranging, 
E r 
~ RPC +L PPC +1 i] 
pc 
pEC 
7 RACs Cat 
PEC 


; cr [+E p+ re] 


ee TT 


i 


(] 


(1) 
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2 


1 R? 
In place of —— write ; this is justifiable if 
L 41% 


* (Goer! 


yes Sant 
( 2s ) Py oe 
It has already been proved that 


p 
pta 


This is comparable with 


[I] = ew! 


p i 
pt (R/2 L) 1] = e-R2L-t 
Thus, 


(p + (R/2L))® 
p E 4. e-R2L-t 


& 


L (6 + (R2L)2 iL 
It will be seen from the foregoing examples that solutions by 
Heaviside’s calculus are simpler and more straightforward than by 


the ordinary calculus. 

It may sometimes 3 
pressions for various parts of the circuit. 
adopted will be illustrated below. 


be necessary to evaluate the voltage ex- 
The method to be 


50 AN INTRODUCTION TO ENGINEERING MATHEMATICS 


With reference to case ) of the example ust discussed : 
- RL. 3 Te 
a= Ee"! —— sinh (Vea R? Be t 
5 re 
L Rie ak R 412 
A ie 


The voltage across the resistance = Ri = ep 


-R/2L.t a ae 
= Re sinh ( / B ae 
Ti ei oh HERES IG 

Aids 7 Tie 


The voltage across the inductance = L di/dt = 


ee ey (ety, 
pape a7 Le oo a ) 


41? 


ie af Ro 
at 412 LC | 
| 
ie iL | 
cosh (Vv ) 
4L? JAE 
Io I 918 PSB ES ; R? 1 
sinh ( = Stic) t | 
Pe) aan 412 LC | 
402° Le 
= E eR cosh ( Seay =) a R.E ehh! Z 
4L? 1G 2 
2L R ipvetlies 
412? LC 
R 
sinh ( piss ) i : 
ee Le 
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Te, is 

= Bese cosh ( — _ —) ee 
41? LC 


: RS il 
sinh ( — -_——)t | : 
41° Le) 
The voltage across the condenser may be found by evaluating ~ 
the expression 


i 
alle: (sees) 
or by re-arranging the voltage equation. 


di , Die » 
ee R = Jdt =E 
+ Pere || 


di 
di : ih ff 
soi eae arae 
E-LO-Ri =— | i 
thus : aa 
cr 2 -R/2L.t 2 1 
4 [ia au - Se sion ( ~~ -=)! 
@ ves 1 mie > Te 
FREE > AEC 
q 2 
- penet+[ can (VE - =z) = Bae ED) 
4L2 LC R2 1 
iV 
AT ~ UG 


N Te Sil 
sinh ( === ==) t 
AVE ILC, 


eft Re-Rek-t va ( 4 EG ; 
— = sin San ere 
eek 412 ia) 


412 LC 
—-— 
— e- RPL. +t cosh (ea) 3 
41? ie, 
Re-R.t SS 
in (VE 
— VOILA SS 4L? IEG, ) 


42) Le 
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R 2 
= e{I — eRMLt Tat aan ca ( R =e aa t 
| it pee a ees IG. 


2 
sinh ( aoe) || 
41? LC 


) In case (2) the voltage expressions are found in exactly the same 
| manner. They are: 


| Zs 9 GEXe 
Ren) as 
| + cosh ( ->)! 
) 41? LC 
R 2 ic 
R 1 sin (YR 1) , 
EH ary. £ 4L? LG ) 
41? LC 
a ie R 
. =E {1-¢Rew. [ cosh (v- ae = met a 
( 4L2 LC Ly ayatlhs 
] 41? 1G; 
| moi 5, 
(sinh (aE sid ( A eee ‘)]} 
412 ae) as are eG) j 
a eae R 
= Eft ener oe - eee ee NT Rag ca 
417. LC ay/ Ces 
4L? LC 


Voltage across resistance = Ri = e, 


R Ee-R2L.t 
= SS sin (vz a mage, 
F Ay ae R2 4L2 


We wane 


Voltage across inductance = L di/dt = e, 


vo R 
= BeR?2L.t [ cos ( +=) Es coe aL are i Re 


LC 412 


a We-= = aa ] 
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1 
Voltage across the condenser = rom fidt=e 
R 


1 R 
=E ‘ ~ ene-t[cos (VE =———)t+” Tq pe 
( We 4S MEAL ARS 


sin — - 
: LG 
Similarly for case (3) 
RE 


£ e-R2L-t 


L 
Ret -R/2L.¢ 
B [1 3 al ; 


eo = {1 - [1 + all e-nntt} 


It is known that the energy stored in an electrostatic field (capa- 
citance) is of the form, W = $C V? 
The corresponding expression for the energy stored in a magnetic 
field (inductance) is, W = +LI? 
possible to evaluate these energy expressions. 


a = 


CS 
Ss 
iH 


Therefore it is 
Table I. summarises in detail the various expressions derived 
from the foregoing discussion on the R.L.C. circuit. 
The time/current voltage characteristics for each case are shown 
in Fig. 12. 
It has already been shown that real equations may have com- 
plex roots. Thus : f 
d (m+jn) =¥ +7), 
where 7 only is unreal. 
Substituting in a typical quadratic type of equation, 
p (m+jn)? +9 (m +jn) tr 
pm? —pn® + Zpjmn + qm +qjn+r 
Pp (m2 —n?) + qm +7 (2pmn + qn) +7 
p (m-jn)yP+ 9 (m -jn) +7 f 
pm? — pn 2Qjpmn +qm—-qyn+r 
P 
$ 


iol 


Similarly, 
(m2 —n®) + qm -j(2pmn +qn) +r 
(m—jn) = x -Jy- 


ONS ae fr | 
[Hf 3] 2a Alf ope R. ay aa ; 
ag ep) al) 


| pis y3 

sah |* i z Hes Kez) Hyde , 
os Oa wat ay 2 ay _2N)| 

; f (3-3) (3 Ei es +4(5-3) (329 ( = 3 ae 


4 233-1 
3° i, Si _ 3}? 712 a -s03] 92 Ha 
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JES 
oak | He} | Sag 


W=).8 
ag? 9783 wed a ah TN] 9 23 wed a | 
“PIst] wosory | ou arerariy) Eacoot 2p 1» % “jooy 


voisupdxy 


pe ee 


SAAS -"919 _ ALINISID 


“ID agVL 


a ee 
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Simple Substitution in the Expansion Theorem, when Roots are 
Complex. 


From A.C. theory it is known that vectors can be represented 
symbolically. The following relationships may be defined : 


Reference to Fig. 13. 


y 


V. Voutace 


AE Current 


Fig. 13.—Symbolic representation of vectors. 
Voltage V may be expressed V/a and current I. may be 
expressed I /p. 


/a means that the vector V is displaced from the X-axis by 
an angle a. 


/@ means that the vector I. is displaced by an angle 8 from 
the X —axis. j 


Vector V may be expressed also : 


V = Va +7 Vv, where V, and V» are referred to rectangular 
cartesian axes OX, OY. and current I = i, -ji, 


They may also have the form : 
Vie =V= Vcosa+jVsina. 
I/p =I =I1cos 8 -jisin§B. 
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Taking a more general case, Fig. 14 illustrates an example 
similar to Fig. 13. 


x + jy may be regarded as being represented by the point Z, 


or as being represented by the rotation of the vector OZ through 
an angle @. 


Fig. 14.—Argand Diagram. 


Referring to polar co-ordinates r 6, 
x+jy =Z=r[cos0+jsin 6] (( >0; -~7<9< 7). | 
The modulus of (x +7y) = 7 
w. (Z| = |x+7y| and @=argZ 
It is known that if y = Ae 


ay = bAe™* = by. 
dx 

Consider the expression, Z = cos 9 + j sin 6. 
& — sin 0 + 7 cos 6. 


= j [cos @ +7 sin 8] 


(ER a 
ay eA 
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By comparison, Z = Aeié 
when 6 = 0 
Z = cos 0 +jsin0 
Be i= il 
| thus I = Aéi® mit Avil, 
ae Z = cos-6+jsin @ = cid 


The following identities exist : 
cos 6+ sin 6 = eid 
r [cos 6 +7 sin 6] = 7 ci0=x + jy 


In general it may be written : 


ze O+2nz eo) 2arc 
aVxt+jy = lavr| [ cos =AB + jsin S42 0] 


a a 
7 (0+2n 7) 
= a vr | € a 
where 7 is an integer. 
If the complex expression for p = -a +78 and the vectorial 


Tepresentation = 7 + jx corresponding to Aei#, then A = V7? 4x2 
and the argument of A = tan-1 r[x. 


Similarly Y = r—jx = A e-io 
We now have : 


A [ei8 el-atip)t + ¢-i9 el-a-ipit] 


3) =) 
( cos f = ao) 
2 
= 2Ae-** cos (pt + 8) 
This expression will now be substituted in the expansion theorem 
in place of 


Yo) 


: b dZe/dp 
where # is expressed as above, 


Cubic Equations. 


Let the general form of the equation be : 
ax? + Bx? + 8 4A =0 (44) 
Reducing this equation in the usual manner (¢.e., substituting 


for x the expression (« <i : <) 
a3 
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If the expression was of the form : ‘ 
x 4+ px + 86x +2A=0 (45) 
the substitution would be (x — §/3). 
With reference to equation (45) the expression becomes : 
B2x Be [ 2 Bx eal 
SB? fae Pa 2 Eo 
[* Bar + oT +B zs + 5 


3 
+o [x-B]+a-0 


= 84% [> -2] + 2 A-35B 42-0 


This is comparable with 


w+ me +n =0 (46) 
(i.e, without a term in x?). 
The Method of solving a Cubic Equation (due to Gardan). 
Let #2) = a0 
i x8 = a3+3a°b +3ab? + 0% 
which, by substitution, may be written : 
a? + 63 + 3abx (47) 
If 38ab=-m (47a) 
and a? +h3 = -1 (47) 


Equations 46 and (47) be similar. 


From equation (47a) : 


G0) = a ore aso = Mas 
3 27 
Consider the quadratic : 
- ms 
2 + nO- 
6 —ni Vw + 4m3/27 
7 2 
Py nu vh ne We 
Gate NiZiiceh Vaya op 


a® and 08 are roots of the equation. 


2 3 4 
r= b= [ page + ae + is ) ] 
x at 2 4 27 


4 ke + (2 . =v 
2) 4 27 
So far it has been seen that for a linear equation 
(ie, ax +b=0; x = ~ dfa) 


there is one root only. 
For a quadratic equation 


(ie., ax®*+ bx +c¢=0, x = 


-b + Ve 4acy 
Ge) 


there are two roots. 


For a cubic equation, therefore, there will be three roots. 


Consider the equation : 
Pail (GG, PS la) 
and let the roots be : 1, a, a. 


The operator j or i (= ~/ —1), has already been discussed. This 
will now be enlarged slightly. Referring to Fig. 15, this gives a 
composite picture of operator j and operator a, the latter being a 
root of the cubic equation under discussion. This has the effect 
of turning a vector through 2 77/3 radians (120°). 


J 
Cane a 


> a°=) 


ie 


Fig. 15.—Illustration of operator’s ‘‘a’’ and ‘‘j:” 
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Writing symbolically, 


aye aes 
Ores el, 2 
2 2, 
= cos Ty j-sin TR @ ei@ms) 
3 
@ = -}3 -j V 3/2 


o 
Noll 


1 (v.e., 2 mradians = 360°). 


Thus, 1 +a+a@=0. 
From the geometry of the figure, it is possible to add two 
intermediate functions of a, namely : ‘ 
at = 60° shift. 
a-* = 300° shift. 


Expressing j in terms of a: : 


= = = 48 

a EL SS: (48) 
: VAS) 

Ba = hei 49 

a ee al (49) 
1 2a 

(48a) 


oes a 
i V3 = Vase 
This operator a is used in the analysis of unsymmetrical three- 
phase circuits. - : 
In th p th  ~ being positi 
n the case 0: Sa te eing positive, 
4 a7 week 
i.e., 27 n2 > 4 m3, there are two complex roots ; and one teal root. 


They are of the form, 


(L) a+b be 
=(240) o. ve 

(II.) 5 De (a—b) 
“-(a+b) °. V3 4 ‘ 

(UL.) 5 I “% (29) ~ 


“In the case of 12/4 + 3/27 being negative, i.e., 27 n® <4m5 
there are three real roots a and 6 have complex conjugate roots 


of the form, 4 ; 
iB 
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» Then’ the roots are : 


(L) -a+v38 
(IL) -a-V3.p 
IIL.) 2a 
a. : n® ms : 
Lastly is the case of equal roots, 7.¢., + =0, all 
1 4 27 
three roots are real and at least two are equal. 
They are: (I) Qa . 
, (IL): -a : 
(IIL.) -& 
In the second case considered, 
n® m3 
——= + < 0 
4 27 


the three real roots yielded, require the extraction of cube roots of 

. complex quantities for a numerical solution. Thus, in its alge- 
braic form, the expression is irreducible, therefore trigonometrical 
solution is required. F 


The Quartic Equation. 


As the degree of an equation becomes Wiener more specialised 
solutions are required. The highest degree of equation which will 
be discussed in this work is that of the fourth power in x, namely 
the quartic or biquadratic equation. 


Method of Solving a Quartic Equation. 


A quartic equation may be of the form : 
m+axe+bx?+cex+d = 0 
This may be reduced to a pair of quadratic equations by finding . 
a real root, mee of the cubic equation : 
— 4 by®+ 2 (ac—4d)y — [c?+ d (a?-4b)] =0 
a four roots may be extracted, by solving the two following 
quadratic equations : : 


xe [4 = V8 m5 + 2y, -)| x + (M+ Reo =10) (50) 


et Ae -VyP-d)=0 (51) 
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The roots are: 


[EV] 
a ge 2 42 #20] an [out V¥E sea || (52) 


7 1- [5+ 2 ies 2y,—b ,]- 
zi) & ae m, Bailes) 
9. + wer 2y,—5 | -4 [o.- yd f (5¢ 
THE SHIFTING OPERATION. 
From the theory of the calculus it can be proved that : 4 
DetA (A=¢ (t)) (54)* 
= et DA + Aaet 
= ett (Dial A (55) 
Similarly, D2 ett A 


= APD eee mead 
= ett [D + a] [D + a]-A 
Thus in general terms : 
DvettA= 4 (D + alt A 
also, @® (D)" ett A =et@ (D+apA (56) 


Heaviside writes his shifting operation (reference “‘Electro- 
magnetic Theory’’) as : 


erat a) 


P fA e-tt] = ent 
[A e-**] ai 

=ie*t[p= a] A ~ (57) 
It is seen that when the operator # is transferred to the right 


of the exponential function e~“* it (p) becomes (p—a). The reverse 
of this process transforms # into (p + a). 


Comparison of Laplacian Transformations with Heaviside’s Methods. 


.The Laplacian Direct Transformation is written : 


Lf = [erro dt (p>0#0 «>. ) 
. “oO . (f() =0 when ¢ < 0f 
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[erria =90) 
L f(t = 6 (0) 
Inverse Transformation. 


The direct transformation may be inverted, 


= 1 pctico 3 * . 
fO=-Z| Lea C>y | 


L* $ (6) = f 


(Note.—The function undergoing transformation’ is sometimes 


denoted’ fi) N)s 
Transformation of unit function, 


5 $(p) = [. e-Pt f (t) dt 


A : 
2 
oe i ert | == 
. p are: 
_ Thus Heaviside’s operation is the reverse of. the Laplacian * 
operation. 


Let A be a function of #, f (t) 
and 6 also be a function of #, ¢ (t). 


Let the term f (t) [1] ¢ (/) be some point on a wave expressing 
a function at a particular time interval. 


Performing a Laplacian transformation on F(X) 1) (8), we have, 
LF (A) -|1) (8) = foe? F (A) (1]-(8) dt (58) 
Jo ; 


Referring the ‘time interval to a new set of Cartesian ordinates 
O’x, O'y, where it is understood the previous ordinates were Ox, 
Oy. . The new time-interval is 4, = ¢-B. 


Thus equation (58) becomes : 


is evt+p) F (As) [1] (8) dt 


0 


= ePp | e-Pt F (4) di 
o 


‘ 
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By definition, (6) = | ePt f(t) dt 
O . 
“ e-Ptep) F (AK) [1] (84) dh = &?B 6p) 
“o 
The factor e-P£ being the shifting operator. 


METHOD OF TREATING ALTERNATING CURRENTS. 


hare, ; eis — e-ix 5 
By definition, sin « = — > - 
; 27 * 


ee 5 [ex — eix] 


It has been sliown that , 
and z z 1 = ett iy? | 
This —2— yy = es! 1) 
ep Lge mee tale 
: == ie te tte] wy -£ raed tn : 
sie a 
Giniady’.,, ose phe en : 
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cosa Pom + 2m] 
3 =) fia | a ‘ 
eae Bee (60 


A more familiar form of alternating voltage is sin (wt + 4). 


‘sin (wt +) [1] = [sin wt cos f + cos wt sin ¢] [1] 


sn wt = Sana cos d = } [ < + eid | ; 
+ w? 
Miko. : : 
cos wt = a2 Ny sin 6 = 2 [ e ~ e-i¢’ | 
pr + w? ‘ 27 
/ ihe ; pw eg esp 
CT es Gora 


| pe ed - =) i 
Pa acces Coit veng) 8 


= pw (cos ¢). + p° (sin ¢) fj 
0," p + w2 ‘ 


Similarly 


2 (cos 4) + pw (sin-d) .’ 


p? + w* 


cos (wi = g) (4) - 2 (0 


When using the foregoing methods it is sometimes found that 
the denominator of the expansion theorem is in the form of a product. 


+ When this contingency arises, the following simplification may be 
employed : 
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Let Zip) be-a product of the form:  fip)-- Fip) 


The roots of fip) -0, are a. b, c. 


The roots of Fi) ze 0, are a, B, A. 
Z' (a) = f° (a) . e (2) 
. 2 (b) = f" (0) . F (0) 
Z' () =f (0) Fo), ete. 
Since f (a) 
f (2) } =v 
f (0) etc 
Z' (a) = f (a) F’ (a) 
Z’ (B) = f (B) F’ (8) 
Z' (A) = f (A) F’ (A) etc 
Since F (a) } - 
' F (g) = 
F (A) etc. 
Then 
5 Ee ( MOG MONG ey ey 
Zo afy Fa bf Fu. ¢ fe Fe 
Yay ge Y(p) eBt Yay eat ) ils 
+ 


Fa. fa F(a) BA) FB) Nf F Fa) 


= K 


E sin.{ot + el 


Fig. 16.—R.L. Circuit. 


Consider a resistance and inductance in series, across which an 
alternating e.m.f. of the form E sin (w¢ + 0,), is applied. See a 16. 


The operational solution is of the form 

E pw cos + p* sin 6, 
R+L ~ p+ w 
S845) 


pw cos 0, + p? sin 4, 
(R+AL) (p°+ o%)) 


Wo 
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Comparable with fip), Fi), © 


Roots of fp) = 0, are 
fi = - R/L (a=f,=) 


Roots of F@ = 0, are 
2 = +jw (a=p.=) 
bs = — jw (B= 2=P) 
fy = 7+ Pl 
Pf -= pL 
Fo) = p? + w? 
pF = 29° 
Pfo Fo = pL (pf + o) 
» Pho Fy = 26? (7 + pL) 
Substituting for 9, : Py 
Tale 7 4 
PL (f?7+ w?) = - = | 
Pea 
2 x P+ otl? | 


From Fig. 17, wL/R = tan" 6, 


Substituting in Yi) 


2 
cos + ry sin 6, 


LRw.cos 6, + R¢sin® , 


L2 ae 
dy R [Rsin 4, — Lw cos. 64] 
= - 
= = [R.sin 6,- Lw.cos 4] 
R 
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Fig. 17.—Pertaining to R.L. circuit. 


aL 


R = VR?+o?L? cos 05 
@L = V R?+?L? sin 6, 
R a 


= [VR?+ w?L? cos 6, sin 4 


— v R?2- w?L?. sin 0, cos 6] 


R 
= 7a [(V Rte? LT) sin (0, - 6.) ] 


The first part of the simplified expansion theorem may be 
written : 


R/L? [ (WR? + @ 12) sin (8, — 83) ] 
— R/L? [R? + w? L?] 
E.e-RiL-tsin (0, -05 


E . e7Ri-t 


The next root can be dealt with in a similar manner. The 
solution is well-known and is of the form : 


fe ana ee a OL 
V Ret L (61) 


To some readers this process may seem rather involved, although 
actually each step is straightforward. However, for the benefit 
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of those readers, the same circuit will be solved by the ordinary 
calculus in an abbreviated form. 


Firstly, the voltage equation is established, and is of the form : 
ines Oe Sey. sin-coe 
dt 


The complementary function gives the transient part of the 
solution, and is : 


(R+LD)i =0 


i = Ae-RILt 

1 di F 

aR ete RE cee 

L a; r E sin wt 

1 d 1 

Gert aint E sin wt 
u E sin wt 
iy wi 

Oo = (D = d/dt) 

= ie Bap 
TRE + 


sin wt may be written in the form: i¥t 
B/Leein yb /Lxext 


“ R/L+D ~~ R/L+ jw 


Rationalising and discarding imaginaries, 


t 


Ree, sin wl —’ —— cos wt 
I ee, 
7 RL? +? 
R. 2 
— sin wt — cos wt 
L2 
i R2 + w?L?2 
L2 


RE sin wt—Lw E cos wt E (R sin wt —Lw cos wt) 
R? + w?L? Fi R? + w? L? 


Pe 
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E R sin wt — Lw cos wt 
VR+e © VR? + wi L? 


R. 


Fig. 18.—Pertaining to R.L. circuit (alternative method). 


From Fig. 18: 


R 
COSs¢) SS | 
V R24 071? | 
"i woL 
cin = | 
V R?+ w? L? 


E cos ¢ sin wt — sin ¢ cos wt 


VR aL? 


E 
eee SIN (OES 
RE GRTR te 
= I sin (wt - ¢) 


Therefore the complete solution is : 
i = I sin (wt-¢) + A.e7RLt 


The value of the constant A depends upon the amplitude of the 
voltage wave at the instant of switching-in. 


When this condition is known the value of A may be substituted 
in the equation, which will then be of a form comparable with 
equation (61). 


DUHAMEL INTEGRAL. 


At a time ?, let an e.m.f. of the form ¢, be impressed upon a 
circuit; shortly afterwards another em.f., et: is impressed at 
time ¢"; after that another e.m.f., én is impressed at time 7, 
and so on. 
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If the operational solution for the first impressed e.m.f. be : f(#), 
then, current due to en = ey f (#). 
current due to én = etn f (&-?"), etc. 
Thus : : 
i= euf (tC) + en fU—0) + con f(O=0") + ew f (-0), etc. 


Let the increase in e.m.f. be small and equal to 8 e, and the time 
interval to $4, 


Then 8i = 2, 8%, f (¢-z) 


If at time, /=0, e.m.f. =eo, then at this time, 
t 


current = é f(t) + | eu f (t-t,) dt, 


° 


If the extreme member of the equation be integrated (hy parts) 
the current solution reduces to a form : 
at 


i=ef(0) + | eh, ae thee dt, 


AL) 
This theorem is sometimes known as the superposition theorem. 


In conclusion, two examples ‘will be discussed to illustrate the 
use of the ordinary calculus and the Heaviside’s operational calculus. 


1 (a). The R.C. circuit (see Fig. 19) 


Ate Sen 
R = Ch =e Ee 
Uae | t 
Differentiating : 
di 1 
= = = O 
a oe 
di 
RiGk=— = 
Cc me 
RC di + dt i = 
i di 
a Rey) 5 
t RC 
Integrating : 


log i + aia + K = O (K & integration constant). 


a = e-HRC | gk = C g-t/RC 


(C is a new constant equal to the constant e*). 
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= i 4b ' 
> \opaeret ere 
ay aac 
tall PVN 


Fig. 19.—R.C. circuit. 


At the instant of switching, there is no voltage across the con- 
denser, so that the current is equal to E/R. 
i = E/R = CetRC = Cee =C 
Thus 7 = E/R e-“RC 


(1b). A similar voltage equation may be written : 


Rao a /G aljP ss. =) 


1 
t ‘= (E [1] 
R + 1/C. 1/P 
Cc 
=E : [1] 
RC#+1 
Therefore : 

Yip) = Cp 

MG = © 

Zo) = RCA +1 

Zo) =p ll 

Zq = RC 


To extract the roots, equate Zip) to zero. 
Ast Ss ROP SO 
pe 1/RC 
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Substituting in the expansion theorem : 
0 Cp ert 
i 1p lee + 2 


pRC 
C (=1/RC) e-WRe-t 
re [o + (-1/RC) RC | 


— 1/R e-URC-t 
as [ 1 ] 


5 
= E e7t/RC 


(2a). The R.L. Circuit. See Fig. 20. 


Fig. 20.—R.L. circuit. 


The voltage equation is : 


a 
a it 
Ri +L ii E sinw 
Dividing by |"_ and re-arranging : 
di IR 13) en 
CE MC eT Baa 


It will be seen that this is a type of linear equation of the form : 


dy fe 
Pe ean 10) 


x 
or 
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Thus : 
GefRiL-dt = iz . sin wt ef/RL-4& + KK dt 
Integrating and re-arranging : 
E g-RiL.t [ese sin wt + K.e-RE-t dt 


eS eho 


The term fe®/"-+ sin . wt . dt must be integrated by parts. 


J eR/L-* sin . wt. dt is comparable with 
|»  at(- uv — Je ae at ) 


Let w = sin wt 
ms dt 


du 


— =ocowt v 


dt 
r fe 


| eR/L-! sin wt dt = sino & eRL.t — RE eR/L-t w cos wt dt 


= EL ktetsinwt = Ke | eR/L-# cos wt dt 
Let = t B 
e u cos w dv Us 
dt 
du ‘ 
oi = - w sin wt 
ah es L eRIL.t 
R 


[a cos wt dt = cos wt x CSS = = eRIL-t (—wsin wl) dt 


L 3 
2 | eR/L-t sin wt dt 


: L ‘ Lw [: 1. Ayaan 
. R/L .t as R/O sinveit =: ===" "costar eek e- 
Je sin wt dt = = eé ote 
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1b; Z 
joe [ eR/IL-t sin wt dt | 
ie 
L?w 
| = x eRL-t sin wt — eRIL-t cos wt 
f L? w? ; 
R | eRL-t sin wt dt 
Pe i ~ i 2 1m ib 
| eRL- sin wt dt + — | eR/L-t sin wt dt = z eR/L«t sin wt 
L’w 
= eRIL- cos wh 
. L2w2: 
eRL- sin wt dt [3 + =| = = eR/L-t sin wt 
Lew 
ar eR/L-t cos wt 


[ eR/L-t sin wt dt = eR/L-t [z sin wt — 
R 
Lw? 
[di een | 


. R Hy 
| exis sin wt dt = *eR/L-t FE sin wt — w cos wt | 


L?w 
cos wt | 
R2 


which may be reduced to : 


| wo? + RL? 
¢ Thus : 
ee E [R/L.sin wt —w cos wt} ea aeRO 
iE; w? + R/L? 
(obviously the terms e~®/L-! and etR/L-! cancel). 
d 193; ne 
a + Ke-RLet 
L (w? + R?/L?) if 
E 


SSS R sin wt — wl cos wt | + Ke-RiL.# 
L? w* + R® 


R sin wt — wl cos wt 


E 
9 VRa Let [ VRP Law? 


| + ke 
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, From previous examples, it may now ie written : 
Pee 
“VR Lat 
where 6 = arc tan wL/R.’ 


gin (wt — 6) + Ke-RiL-t 


The value of the constant K depends upon the position of the 
voltage wavé at the instant of switching-in. Suppose this to be 
at a time T. 


If ¢ is assumed to be zero, when 7=0 and 2=0, 
F di A 
= peSy= t 
[é Ri + ai E sin wi | 


then with time from ¢=0 to 4=T, at ¢=T. 


OO) Sb Tsint (Coy — 0) ee Eee 
Fee TO 
pe =) 

: ( VR + et 
. K=—Isin wT - 6) e+ 
Finally, . 
r | ¢ =I sin (wt-6).- Te“ tt-D sin wT — 6) 


(It is interesting to note that if the instant of switching-in is, 
arranged so that w T = 9, then the transient term is eliminated). 


(2b). 


e = Ri+Lpi =E sin ot = Ee [1] 
pLli+Ri= Ed" [1] 


, E drt 
+ Sip eRe 
oy EM ee Pies 
pL+R "p24 w? 
Equating Z) = 0, Roots = fs = -RIL 
py a jw 
bs = + Iw 


i= Cye-Rh-t + Cy e-ivt 4 °C, et int 
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From elementary electrical theory, 


E sin (wt -— 0 
Z 


Cy et + Cy etivt = 
Z- VR 
@ = arc tan: w L/R. 


Xo = wb: Zy= (R+pL) (p* +e") | 


1 


In accordance with previous discussions on products appearing 
in the denominator, it is permissible to write : 


Pf Fe = PL (p2 + w%) : | 


R R R24 w?L? 
ee | et at || 2 oR |e ee 
L Go| [ L2 ] 
-wR 
ie 


-Y (p = -RJL) = 
Thus : 
-wR/L al wL 
R/L? (R2+0°L?]. VRieel?  VRe G0? L? 


The complete solution is : 


fe 
5. BERGHE = sin (wt —6) e-RiL-t 


Z 


e [ sin (im ome -Rictsin (ot -0) | 


/ 


which is comparable with the previous-solution. ° ' 7 


CONCLUSION. 


Throughout the pamphlet the’ use of the expansion theorem 
has. been stressed. ‘When one becomes familiar with this branch 
of operational mathematics it is often found easier to segregate 
the various functions of into partial fractions (in a similar manner 
to that adopted in the discussion of the expansion theorem), and 
refer to a table of equivalent operators.. 
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“In the text, one of these equivalent operators was proved, 
namely : - 


. pra 
Others are : iI 
Pp -et 
aren DS . 
p E 
iL) esas ALE stem, 
’ Phas e - fe pu etGs,1eLc: 


As an example of equivalent operatbrs, let the 3 following crane 
be considered : 2 p 


“Let an inductance of constant value and ‘a resistance of 
constant value be short-circuited by a switch.” 


The equation is : ‘ 
Ri+ Lp: =Esin @t+a) [1] 


From the section on treatment of alternating currents, the 
substitution using equivalent operators may be made in one “‘step,”” 
viz. : 

E wp cosa + pf sina 


( R+Ip © p+ w mM etc. 
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APPENDIX I. 
The following list of symbols and definitions may be found 


useful in conjunction with some of the examples in the pamphlet, | 


and also as an aid for further reading. * 


M = Bending moment at section considered. 
il 


through centregof area of cross-section. 


NS 
i 


Bending stress at height y from centre of area of cross- 
section( neutral axis). . - , 
. a Young’s modulus (linear elastic modulus) 
linear stress 
5 linear Strain 
R_ =. Radius of curvature of bean at section considered. 


Resilience of a beam due to bending = 1/2ET f M? dx. 


The Suear Force (SF) at a section of a beam is, the algebraic 
sum of all forces to one side of the section. 


The BENDING MomENT (M) at a section of a,beam is the alge- 
braic sum of all moments to one side of the section. 


e = instantaneous value of e.m.f. (=E sin at). 


G 
i} 
° 
4 


= Maximum value of ‘em.f. 
2 mf angular velocity (mechanical and electrical). 
= Instantaneous value of current. ‘ 


= 
8 

"° 
fe} 


= Maximum value of current. 

= Time in seconds (mechanical and electrical). 
Root mean square (R.M.S.) value of e.m.f. 
= R.MS. value of current. 

= Capacitance in farads. 

= Inductance in henries. 

= Resistance in ohms. 

1/wC (Capacitive) \ evn 

wL (Inductive) 


Boat = & bP 
i] 


Reactance = Xe 
Sait x 


Wo 


Moment of inertia of cross-section pertaining to axis’ 
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Z 


WHAIVaAv V A + Hy eM WOK 


= 
|x| 


[2 or al 


Greek letters in common use : i 


te) 


ge FS eyNrIF +A OX HRIDWA 


= Impedance in ohms 
= VR?4 (Xi -Xo)? = V R?+ (wl - 1/wC)? 
=RijX 


= Admittance in mhos (the word ‘ohms’ reversed). 


= Conductance in mhos = y cos ¢ 

= Susceptance in mhos = y sin ¢ 

= TE = VG +B 5 5 
= Angle of lag of current, 

= plus or minus.’ 


= minus or plus. 


= is not equal. 

= less than. 

= greater than. ‘ 
= less than or greater than, but not equal. 
= equals or is greater than. 

= equals or is less’ than. 

‘is identical to.” 


= approaches. 4 
= absolute value, modulus or signless value, 


(i.e, can be + or —). 


= factorial m (1.2.3:4.5..,. -n)- 


pronounced alpha. 
beta. 
gamma. 
delta. 
epsilon. 
theta. 
kappa. 
lambda. 
mu. 

pi (circles only!). 
rho 

5 sigma. 

. phi. 

iD psi. 
omega 
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DETERMINANT. 


A determinant may be simply expresséd, as an arrangement of 
terms in the form of a square which has a value determined by its: 
development into the factorial quantity of its elements, 7.c., 


8, PsAlie ed VE - My | 
%2 Yo. 22 - s Be 
| %3 Ys 23 = e Bs 
%q Ya 24 we * Ba 
Xe Ix 2x - - px 


is a square array (square formation) of x? elements. 


The principle element of the above determinant is : 
t %1, Vo, 23 - - bx . 
+ A minor of an element (say x,) is the determinant formed of the 
remaining elements when the column and row containing x, are 
suppressed. 
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APPENDIX II. 
Example. x 


_ The following example is taken from “Transients in Electric 
Circuits,” by W. B. Coulthard:(Pitman), page 28. 
Sweep Circuit for Cathode-ray Oscilloscope. 


Introduction. 


In most commercial cathode ray oscilloscopes, a sweep or time 
base circuit is built in, to enable current or voltage v.-time wave- 
forms to be analysed. It is usual to apply, the varying voltage 
to the’ vertical deflection plates and simultaneously sweep the spot 
in a horizontal plane. This horizontal base, which is proportional 
to time, is called the sweep. = 


One method is by use of a reactance. 


It is known that growth of current is in accordance with 


[: - oe | 
ati? 

[= ted] 
21 


As the time is here of the order of microseconds (u secs), it is 


seen that 
; BN! Be ne 
ak G) pve plats 


the sweep deflecting force varies inversely as L. 


1 = 


ayes Ol be 


Le. 


For.a particle of charge e moving at a speed v, the equivalent coil 
=ev. IfB is the flux density set up by the sweep magnet coil, . 
the force on a particle is ev B per unit length. 

But force = mass x acceleration, 


and b= tat 


= 
I 
toe 
a 
ae 
red 
I 
we 
AS 
& 
sls 
w 
Se 
ae 
% 


2 
or ZL -()* (Se) <" 
: ja J 2m 


where / is the width of the screen out. 
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APPENDIX III. 
Comparison of Dimensions of’ Units. 
(For use with appropriate Differential Equations). 
Quantity. iigarreret et ELECTRICAL - Heat | 
(Dynamical) (Electro-physical) [ (Thermo-physical) 
Length 3 1) | (} | 
Time, *(T) a fs) (T) 
Mass, [My] my | [M) 
Force, ... ome. | og a) oo 
Energy (Torque) | [M) (L}* (T]-* (] (L)}* (T)-* | 
Quantity of Heat, : [M] TL]? (ry 
Power, ... (M) (L]* [1-3 ™) ye re | 
Electric Charge, (Q] 
Velocity, (L) (7) | 
Current, 1 (Q) (T)+ i 
Pressure, [My (Ly cry? | 
Momentum, (M] (L) [T}" 
Gravitational i | 
Constant, ©... | [MJ-? [L]? [T]-2 | 
|| Voltage, in toy ce (0) ' 
Resistance, [M] [Q]-* (L}? (T]* 4 
Inductance, ‘| oy [Q}-? IL? 
Capacitance, [M}> [Q]? [L]-2 al 
Temperature, ... i ) 
Thermal Capacity, | (Ly? (ry? (a) 
Thermal 3 iM | 
Conductivaty, Oy (1) cars (6) 
Pest Fe | (M] [T]-* [}-! 
Entropy, 


OM) (Ly? fry? (8)2 
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93-4-; 


Helix Angle and Efficiency of Screws and Worms. 
Approximate Radius of Gyration of Various Sections. 
Helical Spring Graphs (Round Wire) it 

cs 3 » (Round Wire) Connected. 

Bs 8 x (Square Wire) J 
Relative Values of Welds to Rivets. 
Ration of Length/Depth of Girders for Stiffness. 
Graphs for Strength of Rectangular Flat Plates of Uniform Thickness. 
Graphs for Deflection of Rectangular Flat Plates of Uniform Thickness. 
Moment of Resistance of Reinforced Concrete Beams. 
Deflection of Leaf Spring. 
Strength of Leaf Spring. 
Chart showing Relationship of Various Hardness Tests. 
Shaft Horse Power and Proportions of Worm Gears. 
Ring with Uniform Internal Load (Tangential Strain) Connected 
Ring with Uniform Internal Load (Tangential Stress) fi. 
Hub Pressed on to Steel Shaft. (Maximum Tangential Stress at Bore 

of Hub). 
Hub Pressed on to Steel Shaft. (Radial Gripping Pressure between 
Hub and Shaft). 

Rotating Disc (Steel) Tangential Bales Connected: 


Ring with Uniform External Load, Tangential Strain \ Connected 


» ” ” a ” Stress 
Viscosity Temperature Chart for Converting Commercial 

to Absolute Viscosities, 5 Connected 
Journal Friction on Bearings, : 
Ring Oil Bearings, 4 
Shearing and Bearing Values for High Tensile Structural 

Steel Shop Rivets, in accordance with B.S.S. No. 

548/1934, Connected. 


Permissible Compressive Stresses for High Tensile Struc- 
tural Steel, manufactured in accordance with B.S.S. 
548/1934. 

Velocity of Flow in Pipes for a Given Delivery 

Delivery of Water in Pipes for a Given Head Connected: 

Working Loads in Mild Steel Pillar Shafts. 

Involute Toothed Gearing Chart. . 

Steam Pipe Design. Chart showing Flow of Steam through Pipes. 

Variation of Suction Lift and Temperature for Centrifugal Pumps. 

Nomograph for Uniformily Distributed Loads on British | 
Standard Beams, | 


Connected. 


Notes on’ Beam Design and on Use of Data Sheets, Nos. 


[ 
84-5-6. | 
Curve Relating Natural Frequency and Deflection, nee 
Vibration Transmissibility Curve for Elastic Suspension 
Instructions and Examples in the Use of Data Sheets, \ 
Nos. 89 and 90. 
Pressure on Sides of Bunker. 
‘5-6-7. Rolled Steel Sections. 


Connected. 


98-99-100. Boiler Safety Valves. 


Orders for Pamphelts and Data Sheets to be sent to the Editor, The 
Draughtsman, cheques and orders being crossed ‘‘Draughtsman 


Publishing Co., Ltd.” 


(Data Sheets are 2d to Members, 4d to others, post free). 


